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I. INTRODUCTION 



has not reached the escape window up to time t obeys 



A generic problem in cellular biochemistry is to esti- 
mate the time - the so-called Narrow Escape Time (NET) 
- that a randomly moving particle spends in a bounded 
domain before it escapes through a small window in the 
domain's boundary. A particle can be an ion, a ligand, 
a molecule, a protein, etc. A confining domain can be 
a cell, a microvesicle, a compartment, an endosome, a 
caveola, a dendritic spine, etc. A variety of processes 
in which the importance of the NET problem is striking 
were discussed in [H-Q- 

Conventional analytical calculations of the NET rely 
on the assumption that the confining surface is perfectly 
reflecting everywhere, except for the escape window - an 
aperture of typical size a. For Brownian motion, evalu- 
ation of the NET probability density function (PDF) Ft 
amounts to finding the solution of the diffusion equation 
with mixed Dirichlet-Neumann boundary conditions jsj. 

In three-dimensions (3D) one finds p] that at suffi- 
ciently large times the probability St that the particle 
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where the symbol "~" signifies that one deals with the 
leading in time asymptotic behavior and omits the nu- 
merical prefactors. The characteristic decay time t^u 
(the subscript "3D" specifies that the search for the es- 
cape window proceeds via the bulk diffusion) in Ecj^. [T] is 
given by 
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where V is the volume of the domain and Dq is the 
bulk diffusion coefficient. This result holds for any 3D 
bounded domain, provided that the boundary is suffi- 
ciently smooth and the ratio a/ R, where R is the typical 
size of the domain, is sufficiently small. 

The PDF Ft than follows via the relation Ft = 
—dSt/dt. Hence, t^D in Eq. [2]can be interpreted as the 
mean time of the first passage to the escape window - 
the mean narrow-escape time. Note that the results in 
Eqs. [THU have been obtained earlier in 

M and i m 

the special case of a sphere of radius R and the escape 
window being a geodesic disc of radius a, a/ R <^ 1. 

To get some idea of typical NET scales, consider an 
example mentioned in [l| - search for the tubule entrance 
in a vesicle by a diffusive ligand. The vesicle size R and 
the radius of the tubule entrance a are of order of 10~^cm 
and 10~^cm, respectively, while the ligand diffusion coef- 
ficient Dq is in the range 10~^ — 10~^cm^/s. Thus t^D is 
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of order of 10~^ to 10~^ seconds, depending on the value 
of the bulk diffusion coefficient. Of course, one may en- 
counter considerably larger first passage times for larger 
R or smaller a, as well as under conditions of molecu- 
lar crowding emerging due to complexity of the cellular 
environment. In the latter case, an effectively subdiffu- 
sive motion can emerge Q. On contrary, interactions of 
particles with molecular motors may induce an effective 
biased motion and thus reduce the NET. 

The analysis based on the "perfectly reflecting wall" 
assumption misses an important factor. In realistic sys- 
tems, in addition to the short-range repulsion, there are 
always some attractive interactions between the surface 
of the domain and the diffusive particle. Capitalizing on 
ideas of Adam and Delbriick, put forward for chemore- 
ception [8| (see also the discussion in [6]), one may sup- 
pose that if such interactions are sufficiently strong, the 
actual search for the escape window will be a two-stage 
process, in which the particle will first find the surface of 
the cell and then will move diffusively along the surface 
until it finds the escape window. Consequently, one may 
expect that in this two-stage process the rate at which 
the escape window is found (and correspondingly, the ap- 
parent narrow escape time tapp) will increase (decrease) 
by an amount that depends on the surface diffusion co- 
efficient. 




FIG. 1. A path of a diffusing particle starting at point "x" 
and leaving the sphere through the escape window placed on 
the south pole. Excursions in the inner part of the sphere 
between two consecutive contacts with the surface (Brownian 
excursions) are marlced by a blue color while the excursions 
along the surface - by a green color. 

It would be even more realistic to suppose that in the 
presence of particle-surface attractive interactions, the 
search for the escape window will be an intermittent two- 
stage process [9-13]: a particle approaching the surface 
will reversibly bind to it and, (if the barrier against lat- 
eral diffusion is smaller than the desorption barrier) , dif- 



fuse over the surface for some (random) time T, after 
which it will desorbs back to the inner part of the con- 
fining domain, approach it again at some other point, 
reversibly bind, diffuse, et cetera. Therefore, as depicted 
in Fig. [l] a typical particle trajectory will consist of a se- 
quence of surface diffusion tours followed by excursions in 
the inner part of the domain, i.e., an intermittent combi- 
nation of diffusion in 2 and 3 spatial dimensions. In this 
case, tapp will also acquire a dependence on the mean 
time Ts of residence on the surface within each round 
of surface diffusion. Given that attractive interactions 
are always present, the result in Eq. [5J based on the as- 
sumption of a perfectly reflecting wall, does not provide a 
reliable estimate of the actual mean narrow-escape time 

^app ■ 

In this paper we study analytically, within a mean- 
field approach, the effect of the particle-surface affinity 
on the mean narrow-escape times. In section II we de- 
fine our model. In section III we first rederive the result 
in Eq. [2] in the spirit of the approach discussed by Berg 
and Purcell \d\ in their analysis of diffusive ligands ad- 
sorption from the extracellular space by cell bound recep- 
tors. Next, we consider a special case of very strong at- 
tractive interactions, such that a particle, once it bumps 
on the surface, stays there for good and wanders along 
the surface until it finds the escape window. This limit 
can be thought off as an analogue of the Adam-Delbriick 
two-stage process. Further on, we consider the general 
case of an intermittent two-stage process. We develop a 
mean-field approach which allows us to derive a general 
result for tapp, valid for any Ts- This result represents 
an interpolation formula from which we recover Eq. [5] 
when Ts — )■ and the result of the Adam-Delbriick-type 
approach for — > oo. Finally, we conclude with a dis- 
cussion of the overall effect of the particle-surface affinity 
on the NET. 



II. MODEL 

We focus here on the simplest geometry in which the 
bounded domain is a sphere of radius R. We base our ap- 
proach on the analysis of representative particle's trajec- 
tories, rather than on the solution of the diffusion equa- 
tion with appropriate boundary conditions. Therefore, 
instead of standard settings " point particle vs window of 
radius a" , we switch here to an equivalent formulation 
in which the particle has radius a, while the escape win- 
dow is a point on the surface. Consequently, we stipulate 
that the particle covers an area 7ra^ , when it touches the 
surface. We suppose, as well, that the particle is initially 
placed at a random position on the surface of a sphere of 
radius rg ~ R — a, a = ja, where 7 is a numerical fac- 
tor of order of unity; its precise value will be discussed 
below. 

Diffusion coefficient of the particle in the bulk inside 
the sphere is Dq. When the particle approaches the 
sphere (i.e., the distance between the particle and the 
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surface of the sphere gets smaller than a), it becomes 
weakly adsorbed and starts to diffuse, with diffusion co- 
efhcient Dg, along the surface of the sphere. At every 
(arbitrary small) time step At, the particle updates its 
state: with probability pd it may detach from the surface 
and diffuse away, and with probability 1 — pd it stays 
adsorbed and continues diffusion over the surface. The 
time T of residence on the surface within a single surface 
diffusion tour is a random variable with distribution 



1 / T 
P.(T) = -exp 



where the mean value 
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Our goal is to define, in this general case, the decay of the 
probability St that the particle has not found the escape 
window up to time t, from which we will define the tail 
of the first-passage distribution Ft and hence, the char- 
acteristic mean narrow-escape time tapp, as a function of 
Ts, Dq, Dg, R and a. 



III. NET PROBLEM WITH A 
PARTICLE-SURFACE AFFINITY 

As the first step, we rederive the results in Eqs. [THU 
adapting to the NET problem the approach discussed by 
Berg and Purcell Q within the context of chemorecep- 
tion. Next, we will extend the developed approach over 
the case when a particle has an affinity to the surface and 
may diffuse along the surface. 



NET problem with a perfectly reflecting wall 
revisited. 



Let Pd = I so that the boundary of the sphere is 
perfectly reflecting. A given path of a particle during 
time t can be then viewed as a sequence of N Brown- 
ian excursions - 3D loops connecting the points where 
the particle has touched the surface; of course, is a 
realization-dependent random variable. These excursions 
correspond to the parts of the path marked in blue in 
Fig. [TJ since we focus on the case pd = 1, the "green" 
parts should shrink to single points on the sphere. 

Note that the term " touching the surface" , as well as 
the hypothetical path depicted in Fig.[l]should be viewed 
with an appropriate caution. As a matter of fact, the 
total number of distinct encounters of a point particle 
with the surface during a finite time interval is infinite in 
the limit of continuous diffusion. To avoid this confusing 
behavior, one has to introduce a finite cut-off distance of 
order of a realistic particle radius. 

Further on, not all encounters with the surface can 
be considered as independent tries in search for the es- 
cape window, but only those Brownian excursions whose 



ends on the surface are separated by a distance greater 
than a; shorter Brownian excursions should be removed 
and considered as a single try. This circumstance has 
been discussed in 'g']. Of course, this criterion is rather 
ambiguous and does not define the precise value of the 
cut-off distance. In this regard, we define a Brownian 
excursion as a part of a particle trajectory which starts 
at a distance a = ja away from the surface and ends up 
on the surface without ever crossing it. In doing so, we 
consider 7 as a fitting parameter which will be chosen 
afterwards in order to match the exact result in Eq. [51 

The probability of not hitting the escape window in a 
single random encounter with the surface is 1— 7ra^/47ri?^. 
If the contacts with the surface can be taken as indepen- 
dent tries, we may estimate the probability that a given 
path, starting at a random location on the surface of a 
sphere of radius R~ a, has not found the escape window 

as a product 1 1 — -jj^ 1 . Consequently, the survival 
probability St will be given by 
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where Pt{N) is the probability that the particle 
"touched" the surface exactly A^ times within time in- 
terval t. 

Suppose now that a particle starting at t = at dis- 
tance cr apart of the surface of the sphere touches the 
surface for the first time at i = ri , for the second time at 
t — Ti + T2, and etc. Then, the probability distribution 
Pt{N) can be defined as 



Pt{N)=E. 



(6) 



where the symbol Er{. ■ ■} denotes averaging with respect 
to the distribution of t- variables, while 9{x) is the Heav- 
iside theta- function which is defined as 6{x) = 1 if x > 
and zero otherwise. 

Using the following representation of the rectangular 
function: 



(7) 



C^^ {. . .} being the inverse Laplace transformation with 
respect to the parameter A, we perform averaging over 
the distribution P{t) of independent, identically dis- 
tributed r- variables and find that Pt{N) obeys 
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is the moment-generating function of the r-variables. 
One may readily notice that Pt{N) in Eq.|S]is normalized, 

To evaluate P{t), and hence, consider the follow- 
ing auxiliary problem - the survival of a particle, (whose 
initial location is uniformly distributed on the surface 
of a sphere of radius R — a), which diffuses with diffu- 
sion coefficient Dq within a sphere of radius R whose 
surface is perfectly adsorbing. The Green's function 
GT-jrlro ^ R — a) solution of this problem is given by 

Q: 



2ttR ^ r 
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Integrating over the angular variables and r we find St - 
the probability that such a particle survives until time r, 
from which we get the desired probability density func- 
tion P(t) = —dSr/dr that the first encounter with the 
surface occurred exactly at time moment r: 
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Note that the distribution in Eq. [TT] has been previously 
obtained in fl5| within a different context. 

Before we proceed further, several remarks concerning 
the PDF in Eq.[TT]are to be made. The distribution P{t) 
involves three different time scales. The smallest one cor- 
responds to the most probable value ~ a^/Do, which 
means that most of the time the particle simply bounces 
onto the surface almost immediately without leaving it 
for any considerable distance. Further on, at interme- 
diate scales the distribution P{t) has a "fat" algebraic 
tail P{t) ~ r~'^/^. In this regime P{t) describes the 
probability for a random walk, commencing at a plane 
bounding an infinite 3D system, to return back to the 
plane for the first time after t steps. As a matter of fact, 
the mean r - the mean length of Brownian excursions Tf, 
- is dominated by this very regime: 



Tb = dr tP{t) 
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and is R/a times larger than the most probable return 
time. Finally, at times of order R^/Dq, finite-size effects 
dominate and the distribution P{t) decays exponentially 
with time. 

Consequently, the moment-generating function of a 
random variable r obeys 
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from which equation we find that at sufficiently large 
times t, the distribution function Pt{N) follows: 
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This distribution is centered around the mean value N = 
t/Tf, and, at fixed i, decays exponentially with N on both 
sides of the N . 

Now, the asymptotic decay form of St in Eq. [S] can 
be determined in two different ways. We can either con- 
vert the sum into an integral and use the asymptotic dis- 
tribution in Eq. 1141 or perform summation exactly and 
then invert the Laplace transform in the asymptotic limit 
t ^ oo (A — 0). We proceed with the latter scenario. 
Plugging Pt{N) given by Eq.[8]into Eq.[5]and performing 
summation over N, we get 



St - c-^ 



(1-0a) 
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In the large-i limit, the integral in Eq. \T5\ is dominated 
by the behavior of (jjx in the vicinity of A = 0. Expanding 
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we find that the asymptotic behavior of St in Eq[Tn] fol- 
lows 



exp 
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Choosing now 7 — 7r/4, we see that the latter decay form 
coincides with the result in Eqs[TH21 Consequently, we 
may interpret t^o in Eq. [2] as 



4i?2 



(18) 



where the first multiplier determines the mean number 
of independent tries necessary to find the location of the 
escape window, while the second factor is the mean time 
separating independent tries - the mean length of a Brow- 
nian excursion. 



B. NET problem with particle-surface affinity: 
Adam-Delbriick-type two-stage process. 

In this subsection we consider an opposite extreme 
case, i.e. that of = 0, so that once a particle happens 
to approach the surface of the domain, it stays there for 
good and wanders along the surface until it finds the es- 
cape window. In a sense, this is an idealized situation. 
Indeed, in this case the barrier against the desorption 
should be infinitely large, and consequently, the barrier 
against the lateral diffusion should be infinitely large too. 
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effectively suppressing the movement of tlie particle along 
the surface. 

Neglecting the time Tb it will take, on average, for the 
particle to arrive for the first time at some random point 
on the surface of the domain, we write the probability 
that a particle diffusing along the surface with diffusion 
coefficient Dg won't find the escape window until time t 
as 



(19) 



where A{t) is the mean area swept on the surface of a 
sphere of radius i? by a diffusive disc of radius a until 
time t - a two-dimensional analog a Wiener sausage. This 
area is defined by a series 
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where 



ak 



dxP,,{x)] I \ dxPl(x\ (21) 



Pvix) being the Legendre functions, while = — 1 + 
tt^ I2B? and Vk are the roots (numbered in the ascending 
order) of the equation 



(22) 



Note that the expansion in Eq. [20] differs by a factor 
(^1 ^ jj^^ from the formal solution of the trapping prob- 
lem on the surface of a sphere flGj. This difference orig- 
inates from different initial conditions. Namely, in our 
case a particle can be initially located at any point on the 
surface (including the area covered by the trap, in which 
case it disappears instantaneously, - this corresponds to 
finding the escape window at a first try), while in the sit- 
uation studied in [l6| the particle starts from a random 
point somewhere outside the trap. 

The leading asymptotic behavior of St in Eq. [iH] is 
dominated by the smallest root vi of Eq. (251 Hence, the 



asymptotic behavior of St is 



St ~ exp 



with 



t2D 
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(24) 



where the subscript "2D" signifies that the search for 
the escape window proceeds in this case via the surface 
diffusion. 

When a/R <C 1, for the smallest root one gets vi 
1/2 In (2i?/a) and, consequently, t2D obeys, in the leading 
order in a 



ho 



2i?^ 



(25) 



This result has been also obtained in [l[ and @, and 
earlier by Bloomfield and Prager [l3| in their calculation 
of the attachment rate of tail fibers to bacteriophages. 
Note that t2D ^ ln(i?) and thus should be much larger 
than Tf, ^ R, which describes the mean time necessary 
to reach the surface of the domain. This means that it 
was quite legitimate to discard this contribution in our 
analysis. 



C. NET problem with particle-surface afflnity: an 
intermittent two-stage process. 

We turn finally to the general case when the detach- 
ment probability < < 1, so that a particle, when 
touching the surface, will remain weakly adsorbed and 
wander on the surface for some random time T, then de- 
tach and diffuse in the bulk, re-attach to the surface, and 
etc. 

Consider a path starting at a random point on the sur- 
face of a sphere of radius R — ira/A and suppose that this 
path touched the surface of the sphere at time moment 
Ti, then wandered along the surface for a random time 
Ti , detached from the surface at time moment ti + Ti , 
subsequently returned to the surface at time moment 
Ti + Ti -I- T2, etc. Then, assuming that subsequent vis- 
its of the surface can be considered as independent tries 
in search for the exit, the probability St that the escape 
window has not been found by such a path comprising 
N rounds of Brownian excursions followed by subsequent 
surface diffusion tours, can be written as 



5. = p.(A^ = o)-HEMn(i-|§)M*-E-^ + E^^0M-^+^^^ + ^^0}' 

N=l k=l ^ ^ \ k=l k=l / \ k=l k=l / 

I 



where now the symbol Et^t{- ■ •} denotes now the averag- ing with respect to both the distribution of r-variables. 
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Eq. [TTl and the distribution of T-variables, Eq. [3l 

Note that Eql^B] taeitly assumes that each surface dif- 
fusion tour is an independent try in search for the escape 
window, which manifests itself in the decouphng of the 
average of the product into the product of average values. 
This is, of course, an uncontrollable assumption. 

On one hand, the distribution P{d) of the distance d 
between the point where the particle detaches from the 
surface and the point where it re-attaches to the surface 
again after an excursion in the bulk is given by the Pois- 
son kernel for a three-dimensional ball, 18]: P{d) ^ 
This is a broad distribution, such that the areas visited 
on the surface in two consecutive surface diffusion tours 
will not, on average, significantly overlap. 

On the other hand, the surface of the domain is of a 
finite extent and one will certainly have an oversampling 
- some parts of the surface will be visited many times 
before the escape window is found. This will incur some 
correlations in the search process, since the true survival 
probability St accounts only for the actual area swept on 
the surface by a particle up to time t, and counts multiple 
visits to the same place as a single try. In this sense, Eq. 
^Sidelines a lower bound on the true survival probability, 
precisely in the same way as the Rosenstock (or Smolu- 
chowski) approximation defines a rigorous lower bound 
on the decay function for the trapping problem (see, e.g., 
0). 

Consequently, decoupling of correlations defines a rig- 
orous lower bound on the mean narrow-escape time. 
Given that, as we proceed to show, tapp obtained via such 
a mean-field approach entails exact results for Ts — > and 
Ts — > c», one may judge that it is a useful and plausible 
approximation. A further discussion of this matter goes 
beyond the scope of the current paper and will be ad- 
dressed both analytically and numerically elsewhere [13] . 

Using the Laplace transform representation of the rect- 
angular function, Eq. [71 we may conveniently rewrite 
Eq. and perform straightforwardly averaging over the 
distributions of r and T-variables. In doing so, we get 



St=C- 



A 



where 



(27) 



(28) 



Explicitly, F\ is given by a series 



Fx={l- 



fc=i 



^afe (l + i^kii^k + 1) 
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(29) 

Consider now the behavior of the coefficients in this series 
in more detail. First, for small a/ R a good approximate 
solution of Eq. [22] is i^fc = fc - 1 + l/21n(2i?/a) [11], 
and thus the roots of Eq. [22] grow linearly with k. Sec- 
ond, in the leading in a/R order, ak — where Sk,i 
is the Kronecker-delta (ak = 1 for fc = 1 and zero oth- 
erwise) , correction terms to this dependence are of order 



of {a/R^fk, where fk is a rapidly decaying function of 
k. All this permits us, so far as we are interested in the 
small-A (large-i) limit, to consider only the first term in 
the series in Eg. 1291 and skip the remaining terms, giving 



Fx 



l + vi{vi + 1) 



DsTs 
i?2 



(30) 

Finally, notice that the first term in Eqs. [26l [27] decays 
rapidly in the time t domain, compared to the second one 
(indeed, its characteristic decay time is just Tf,) and thus 
its contribution is negligible in the large-i limit, we get 

^ ^_^\ ( iD^ia^ +^vi{vi + l)DsTs) , , 
\y 12DoR^Ts+naR3 + 

'(31) 

which describes the asymptotic behavior of the survival 
probability in the limit A -> 0. 



IV. RESULTS AND DISCUSSION 



Inverting Eq. 1311 we obtain our main results: 



St ~ exp 



where 



UDnR^^T, 



(32) 



(33) 



is the mean NET for diffusion in a sphere with particle- 
surface affinity. Note that for Tg —?' Owe recover the exact 
result in Eq. [21 while for Ts we find from Eq. [33] the 
exact result in Eq. [24l specific for the Adam-Delbriick- 
type two-stage search process. We recall, as well, that 
Eq. 1331 defines a rigorous lower bound on the actual mean 
narrow-escape time for any Ts ■ 

It is expedient to cast the result in Eq. [33] into a phys- 
ically meaningful form, using the characteristic times in- 
troduced in Eqs. H [H [22 



n 



^app 



(34) 
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t2D 



This equation has a transparent physical meaning: the 
numerator on the right-hand-side of Eq. [M] defines the 
overall time spent, on average, in a single Brownian ex- 
cursion in the bulk followed by a surface diffusion tour, 
while the denominator defines the average fraction of the 
sphere surface covered within a single tour of surface dif- 
fusion. Therefore, tapp equals the time consumed by a 
Brownian excursion and a single surface diffusion tour, 
times the number of tries necessary to cover the whole 
surface. We note that a similar result has been obtained 
in [2l|, ll^l for search by a diffusive protein for a specific 
binding site on a DNA molecule. 
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Further on, we rewrite Eq. |34] formally as 



t. 



app 



hD[i + -]/[i + 



(35) 



One notices next that tapp is a nionotonically increasing 
(decreasing) function of if t2D > tsD {ho < Hd)- 
Consequently, tapp > tsD when 



Do TT 



R 



(36) 



Eq. [36] is of order of unity and the inequality in Eq. [36] 
evidently fulfills. This means that for this example the 
particle-surface affinity will generally lead to larger mean 
narrow-escape times, compared to the estimate based on 
the assumption of a perfectly refiecting wall. To inverse 
the inequality, one will need the ratio R/a to be or order 
of 10'^ to 10^, which may be realized, say, for catalytic 
reactions in microporous media {R being the radius of 
a pore and a - radius of a catalytic site). In this case, 
indeed, one may expect that particle-surface affinity will 
reduce the effective times of the first-passage to the cat- 
alytic site and thus enhance the reaction rate. 



Typically, Ds is less by two or three orders of magnitude 
than Dq (may be even less under the conditions of molec- 
ular crowding [23]), which means that the ratio on the 
left-hand-side of the inequality in Eq.[36]is of order of 10^ 
to 10^. For the example, mentioned in the Introduction, 
- search for the tubule entrance in a vesicle by a diffusive 
ligand, one has R/a ^ 10^ so that the right-hand-side of 
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